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I.  INTRODUCTION 


Lamb's  [l]  half  space  problems  have  been  the  subject  of  many 
investigations  because  they  are  fundamental  to  numerous  two-dimensional  ap- 
plications. Kraut  [2]  extended  Lamb's  solution  for  a line  source  on  an  iso- 
tropic half  space  to  a transversely  isotropic  half  space  by  developing  a solu- 
tion with  the  Cagniard-deHoop  technique.  A summary  of  this  solution  and  a 
list  of  related  references  are  available  in  Kraut's  review  article  [3].  Kraut's 
solution  is  the  nondispersive  solution  to  the  present  problem  corresponding  to 
a replacement  of  the  laminated  composite  with  a homogeneous  transversely 
isotropic  material. 

In  order  to  construct  a solution  that  includes  the  effect  of  geometrical 
dispersion,  we  begin  with  the  equations  of  motion  and  boundary  conditions 
from  the  effective  stiffness  theory  advanced  by  Sun,  Achenbach,  and 
Herrmann  [4].  One  of  several  other  approximate  theories  that  are  available 
for  composites  could  have  been  selected  for  this  study;  however,  the  form  of 
the  fundamental  solutions  obtained  here  would  not  be  changed  by  such  a choice. 
Two  related  problems  that  demonstrated  the  effect  of  geometrical  dispersion 
on  the  dynamic  response  of  a laminated  composite  have  recently  been  solved 
[5,  6]  within  the  context  of  this  effective  stiffness  theory,  and  those  solutions 
provide  the  background  for  this  investigation.  The  first  problem  [5]  was  con- 
cerned with  one-dimensional  pulse  propagation  in  an  obliquely  Icuninatcd  com- 
posite: the  second  [6]  resulted  in  a steady-  tatc  solution  for  loads  moving  over 
a composite  half  space.  In  each  case,  a far-field  solution  was  developed  that 
contained  low-frequency  long  waves  consistent  with  the  basic  assumptions  in 
the  effective  stiffness  theory.  A near-field  solution  would  require  high- 
frequency  short  wave  information  that  is  available  only  from  the  exact  theory. 
In  other  words,  the  approximations  that  were  used  to  construct  the  solutions 
matched  the  approximations  inherent  in  the  field  equations  and  boundary 
conditions  of  the  effective  stiffness  theory.  The  same  approach  is  followed 
in  the  present  work. 
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In  the  following  analysis,  Fourier  and  Laplace  transforms  are  used  to 
obtain  a formal  solution  to  the  problem.  A far-field  solution  is  developed 
for  the  body  wave  responses  in  the  interior  of  the  half  space  by  first  using 
the  method  of  steepest  descents  to  approximate  the  inversion  integrals  near 
the  appropriate  saddle  points.  Then,  the  phase  velocities  are  expanded  for 
low  frequency,  which  permits  the  inversion  of  the  Laplace  transforms.  The 
head  wave  and  Rayleigh  wave  are  treated  in  a similar  manner. 

The  final  results  show  that  from  two  to  five  separate  waves  contribute 
to  the  response,  depending  on  the  location  of  the  observer  in  the  laminated 
half  space.  The  quasi-  dilatational  wave  and  head  wave  are  shown  to  be  one- 
sided waves,  whereas  the  quasi-equivoluminal  wave  is  one-sided  until  a 
critical  angle  is  exceeded,  and  then  it  becomes  two-sided.  Also,  cusps  may 
develop  on  the  quasi-equivoluminal  wave.  Geometric  dispersion  caused  by 
the  laminations  is  shown  to  strongly  influence  the  body  wave  results  by  elimi- 
nating the  singular  behavior  that  exists  for  the  displacements  at  the  wave- 
fronts  in  the  nondispersive  solution.  The  singularity  in  the  Rayleigh  surface 
wave  is  also  found  to  be  eliminated  by  geometric  dispersion. 


II.  STATEMENT  OF  PROBLEM 


The  periodically  leiminated  half  space  subjected  to  the  line  impulse 
PQfi(t)5(5«j)  is  shown  in  Fig.  1.  Pq  is  the  magnitude  of  the  line  impulse;  6(t)  and 
6(Xj)  represent  Dirac  delta  functions  with  arguments  time  t and  spatial  co- 
ordinate Xj.  Because  the  half  space  is  in  a state  of  plane  strain,  the  equations 
of  motion  from  the  effective  stiffness  theory  [4]  can  be  written  as 


■^1^1122^1  ^3^22'^i  ^4^12^2  ^5^112*^21 

■^6®2*^21  ^7®1*^22  ~ "^1^22^1  ^2^1  ^3^2^21 


(1) 


■^8^1122^2  ^9^11^2  ^10^22^2  ^4^12^1  ^11^1^21 


^12^112'^22  ■ ^13°2*^22  " "'^1®22^2  ^ ^2*^2  ^3^2*^22 


(2) 


^5®112^1  ■ ^11^1^2  ■ ^14^11^21  ^15*^21 


(3) 


^12^112^2  “ ^13^2^2  " ^16®11*^22  ^il’^22 


~ ^3®2^2  ■ ^4^*22 


(4) 


where  u^  and  u^  are  the  average  displacements  in  the  x^  and  yi^  directions, 

respectively,  and  il>2i  ^22  microdeformations  in  the  stiff  leiminations , 

The  notation  a.u,  = 3u./ax.  is  used,  and  dots  above  variables  denote  partial 
1 J J 1 

differentiation  with  respect  to  time.  The  constants  in  the  above  equations  of 
motion  are  defined  as  follows: 


= Tl{\^  + -Ml  - Tl)(\^^  ■;■ 

a3=f.n/(l  -r,) 


a = X + fi 
4 m in 


^5  " '^1 


^6  = ^^3 


a,  = 7]{\.  - X ) 

I ' t m 


a^  = d"  a,/12 

8 in  3 


“q  --  ’I'^f  " " - 


",o  = 'V  * - ■^' 


‘‘ii  =’l'*‘r-''m> 


^12  “ ^^8 


aj3  =HajQ 


^14  = d^a2/12 


= Tlfd  - T])t^(.  + T)8j^J/(1  - V) 


a,,  -d^a9/12 


T]f(l  - T1)(X^  + 2h^)  + T1(X^^  + - ’l) 


(5) 

( Cont. ) 


Utoi  ( W'> ' 


=T1P^+  (1  - Tl)p^ 

^3  = T]b^ 


b4  = d^b^/lZ 


where 


T1  = d^/h 


d + d 
f m 


The  preceding  definitions  depend  on  the  Lame  parameters,  thicknesses,  and 

densities  of  the  soft  (m)  and  stiff  (f)  laminations,  which  are  denoted  bv  P 

m m’ 

d^,  and  p^  and  \^,  d^,  and  p^,  respectively. 

The  half  space  is  initially  at  rest,  i.e.,  u^,  u^,  u^,  ...  vanish 

at  t = 0,  and  the  tractions  on  the  surface  are  specified  as 


t^  = 0 


where  t^  and  t^  are  shear  and  normal  tractions,  respectively.  It  is  also  re- 
quired that  all  displacements  and  stresses  vanish  at  an  infinite  distance  from 
the  load. 
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III.  FORMAL  SOLUTION  OF  PROBLEM 


In  order  to  develop  a formal  solution  to  this  problem,  an  exponential 
Fourier  transform  defined  as 


/a> 

f(Xj)  exp(iKXj)dXj 


with  inverse 


f(Xf)  = ^(k)  exp(-i«x^)d« 


and  a Laplace  transform  defined  as 


f(o’)  = / f(t)  exp(-at)dt 


with  inverse 


f(t)  = J f(a)  exp(o-t)dCT 
*^Br , 


are  applied  to  equations  (l)-(4)  subject  to  the  quiescent  initial  conditions. 
Nondimensional  parameters  are  then  introduced,  and  the  resulting  set  of 
equations  is  reduced  to 


B,  - isB  ~ - s^B.n  = 0 
1 dy^  ^ ^ 


2,^  ^ 

B.  - isB,  4^  - s‘^B-v  = 0 
4 3 dy  5 


where  nondimensiona  1 quantities  are  defined  as 


u=Uj/h  , V = u^/h  , x=x^/h  , y = x^/h 


The  coefficients  in  equations  (13)  and  (14)  are  given  by 


where 


Bf  - Aj  - A^/Ajq 


^2  “ "^2  ■ ^5^^11 


B3  - A3  + AgAg/A^j  + A^Ag/A^g 


- ^6  ■ ^9^^11 


B5  - A?  - Ag/A^g 


Ai  = (1  + q\^)(l  - T])s^/12  + 1/(1  - 71) 

= (710  + 1 - 71)  +[ll\6j  + (1  -7l)6^1q^ 


A3  = (*r„+  l)q 


A4  = TIAj 


A5  = TKe^y  - 0^^)q 


(18) 
(Cont.  ) 
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Ag  = Tl(v  - l)q 


Ag  = 'HA^ 

Ajq  = + T)v  + Tl^/(1  - Tl) 

A,j  = s^Tl^A^/12  + 11^6^  + T]^6^/(1  - T]) 


®=Pf/Pm 


6£=€^+2=(\j+  2p.^)/fx^ 


6 = £ + 2 = (\  + 2p.  ) /^ 

mm  m mm. 


Also,  in  anticipation  of  the  form  of  the  final  solution,  5 has  been  eliminated 
by  substitution  of  a nondimensional  slowness  defined  as 

q = f/s  . (2 

Note  that  s is  confined  to  the  dispersive  terms  (those  terms  that  contain  the 
number  12)  in  the  definitions  for  the  A's. 

Equations  (13)  and  (14)  are  combined  to  yield 


®1®4  - <®4®4  * Bi®5  ■ ^ = 0 . 

dy  dy 
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Then,  because  the  general  solution  of  equation  (21)  is  proportional  to 
exp(-msy),  we  obtain 


4>(m,  q,  s)  - m 


B1B4  - 


(^2®4 


B1B5- 


B^)m^ 


= 0. 


(22) 


Since  y > 0 and  s can  be  considered  positive,  we  choose  the  two  roots  of  this 
equation  that  have  positive  real  parts  in  order  to  satisfy  the  conditions  at 
infinity  (y  -•  “) . The  transformed  displacements  and  microdeformations  then 
become 


where 


a = C^e^  + C^e^ 


V = CjDjej  + C^D^e^ 


~ ^2^2®2 

^22  ""  ^2^2®2 


(23) 


e^  = exp(-m^sy)  , j = 1,  2 (24) 

and  Cj  and  are  arbitrary  constants  to  be  determined  from  the  boundary 
conditions.  The  other  coefficients  in  equation  (23)  are  defined  as 

Dj  = (B2  - Bjm?)/im^B2 
Ej  = -m.A4)/AjQ 

Fj  = s(iA5  - >^jDjAg)/Ai,  ' j = 1,  2 . (25) 
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Expressions  for  the  two  tractions  are  available  [7]  and  their  transforms 
reduce  to 

+ ^2«2"2 

where 

G.  = -iqsD.  - m.sA.  - A.E. 

J J J 1 4 j 

H.  = -iqs  e - m.sD.A,  - A„F.  , j = 1,  2 . (28) 


In  order  to  satisfy  the  boundary  conditions  specified  by  equations  (7)  and  (8) 


Cj  = PG2/A 

(29) 

C2  = -PG^/A 

(30) 

in  which 

A = G^H^  - 

(31) 

and  a nondimensional  line  impulse  is  introduced  as 

1 /2  2 

P = P-(Hi  /p  ) ' /p  h^  . 

(32) 

Double  inversion  integrals  for  the  transformed  quantiti 
now  be  written  by  using  equations  (10)  and  (12).  This  ( 

es  in  equation  (23)  can 
completes  the  formal 

solution  to  the  problem. 
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IV.  FAR-FIELD  SOLUTION  FOR  BODY  WAVES 


Since  an  explicit  ev'aluation  of  the  double  inversion  integrals  would 
be  quite  complicated,  if  not  impossible,  approximations  are  sought  in  order 
to  construct  a far-field  solution.  The  Laplace-transformed  displacements 
are  given  by 


We  first  consider  equation  (22),  which  is  a polynomial  with  three 
variables  that  provides  the  roots  m^  and  m^.  If  s is  set  to  zero,  the  re- 
sulting equation  defines  the  nondi sper sive  two-dimensional  slowness  surface 
for  a transversely  isotropic  material  in  plane  strain 


This  slowness  surface  is  shown  in  Fig.  2 for  y = 100,  T)  = 0.  8,  0=3,  = 

0.  35,  and  = 0.  3,  where  the  limiting  values  are  the  reciprocals  of  the  wave 
speeds 


Sf 
m I 


22  j (1  - T1  + TlGHd  - Tl)  + T]6^]  j 


B =6  - ; Y i 

12  21  ~ |(i  . Tl  + T10)[>^(1  - 71)  + T)]| 


For  the  present  data,  3 ^ ^ = 9.  47,  322”  ® 12  ” 

Points  on  the  wave  surface  satisfy 

iqx  + my  - t = 0 (39) 

and  since  a vector  from  the  origin  to  a point  on  the  wave  surface  is  parallel 
to  the  normal  to  the  slowness  surface  (Fig.  3) 


X = -iK$,. 


y = 


where  K is  a constant,  and  the  comma  denotes  partial  differentiation.  Substi- 
tution of  equation  (40)  into  (39)  yields 


K = T/(q$  + m4>  ) 

^ 0, q 0, m 


and  therefore,  from  equation  (40),  the  wave  surface  is  defined  by 


X = -iT$  ) 

0,  q 0,  q 0,  m' 


y = T$  /(q^r>  + rn$.  ) 

^ 0,m''^  0,  q 0,  m' 


In  addition,  m and  q must  be  on  the  slowness  surface,  as  given  by  equation 
(35).  Expressions  for  the  derivatives  that  appear  in  equations  (42)  and  (43) 
are  given  in  Appendix  A;  the  wave  surface  that  corresponds  to  the  slowness 
surface  in  Fig.  2 is  shown  in  Fig.  4. 
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With  the  foregoing  as  background,  we  return  to  the  evaluation  of  the 
integrals  in  equations  (33)  and  (34).  A branch  point  occurs  whenever  m = 0, 
which  reduces  equation  (22)  to 


= 0 


(44) 


and,  for  the  nondisper sive  case. 


2„2 


= (Tie  + 1 - Tl)(l  + q P j j) 


(45) 


= (Tie  + 1 - Tl)(l  + . 


(46) 


Therefore,  equation  (44)  is  satisfied  at  the  four  points 


q = -"  i/p  j ^ 12  ■ 

Other  branch  points  occur  when  these  can  be  determined  by 

setting  the  discriminant  of  equation  (22)  to  zero.  Rayleigh  poles  (q  = T i/3  ) 

also  exist  when  A=  0,  [7].  The  branch  points  and  Rayleigh  pole  located  in 
the  lower  half  of  the  q-plane  are  shown  in  Fig.  5.  Branch  cuts  are  also  indi- 
cated in  the  figure.  The  path  of  integration  along  the  real  q axis  can  be  re- 
placed with  a Cagniard-deHoop  path  in  the  lower  half  of  the  q-plane  for  x > 0. 
This  problem  was  solved  for  the  nondisper  sive  case  by  Kraut  [2]  and  provides 
the  first  approximation  to  the  present  dispersive  probl  jm.  Kraut's  solution 
is  exact  and  valid  everywhere,  including  near  the  wavefronts  shown  in  Fig.  4. 

He  also  determined  that  the  displacements  were  singular  at  the  wavefronts 

-1/2 

with  strength  0(  i t - T=l<i  ),  whereas,  at  the  cusp  tips,  the  singularity 

-2/3 

became  0(  | t - | ),  where  corresponds  to  the  arrival  time  of  the 

wavefront.  These  singularities  are  directly  related  to  the  order  of  the 
- corresponding  saddle  point.  In  the  present  analysis,  the  low-frequency 
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behavior  of  these  saddle  points  is  used  to  develop  a dispersive  solution  near 
the  wavefronts  that  exist  in  the  nondisper sive  solution.  Therefore,  the  fol- 
lowing developnaent  is  valid  only  for  large  distances  from  the  load,  where  the 
response  is  expected  to  contain  only  low-frequency  components.  For  the 
method  of  steepest  descents,  it  was  convenient  to  refer  to  Miklowitz  [8],  who 
presented  analogous  approximations  for  regular  and  irregular  wavefront  cases 
for  a homogeneous  half  space. 

In  order  to  proceed,  we  substitute 

X = r sin  (p  , y = r cos  0 , 0 < 0 < tr/Z  (48) 

into  the  first  integral  of  equation  (33)  and  obtain 

I = ^ exp(-rsf^)dq  (49) 

-CO 


whe  re 


f,  = iq  sin  0 + m.  cos  0 


(50) 


The  method  of  steepest  descents  requires  that 


afj  8m^ 

-- — = i sin  0 + — - — cos  0=0 

aq  ^ aq  ^ 


:51) 


and  assuming  that  cos  0^0,  which  excludes  the  surface  response,  we  obtain 


am^ 

9q 


+ i tan  0=0. 


(52) 
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A value  of  q that  satisfies  this  equation  is  a saddle  point  qj^,  and  the  integral 
can  be  approximated  in  the  vicinity  of  this  point  to  yield 


'{sG^A- 


where 

2 2 

a f,  d m 

f " = 5-  = ^ cos  ^ 0 . (54) 

9q  3q 

This  approximation  to  the  integral  is  valid  as  long  as  the  saddle  point  is  not 

2 2 

on  a branch  cut.  The  calculational  procedure  to  obtain  8 m/9q  is  given  in 
Appendix  B . 

To  invert  into  the  time  domain,  we  first  substitute  (see  Fig.  3) 

m^  = 3 ^ cos  cCj  , *51“  ^ (55) 


into  equation  (50)  to  obtain  from  equation  (53) 

1 (2Trrs)~  ^ ^^(  sG^A  S ! f I ^ exp[-  s(x  sinQ'j  + ycoso'j)/3].  (56) 

We  then  expand  the  slowness  in  powers  of  s (Appendix  C) 

1/3  = l/3j  + CiS^/3j  (57) 
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and  use  the  convolution  theorem  to  invert  the  transform.  The  final  results 


where 


u(x,  y,  t)/P  = " :-fj  Re  lim 

s-0 


L 11'  q=qj 

X ( p-^/^Ai(±6  )dp  - lim(-^\ 

^0  s-^ol  ^ /X2lf"i^ 

X j p“  ^'^^Ai(i5  2)dp 

r 

= 3iC-'(x  sin  o.  + y cos  oJ/P. 

J L J J J J j 

3.  = T-p-(x  sin  O'.  + y cos  oJ/P.  /X- 

<3  L 3 3 3j  3 


in  which  j = 1,  2,  and  the  positive  sign  in  the  argument  of  the  Airy  function  is 
used  if  C-  is  positive;  otherwise  the  negative  sign  is  used.  Expressions  for 
the  V displacement  are  obtained  by  inserting  in  the  first  term  of  equation 
(58)  and  in  the  second  according  to  equation  (34). 

Figure  6 shows  the  locations  of  the  saddle  points  that  correspond  to  the 
slowness  surface  in  Fig.  2.  The  quasi-dilatational  branch  (m^)  has  only  one 
saddle  point  and  it  is  always  regular  because  it  does  not  fall  on  a branch  cut 
for  0 < <p  < tt/2.  If  0 < 0 < 0^,  where  0^  is  the  value  of  0 at  point  2 in  Fig.  6, 
the  quasi-equivoluminal  branch  (m^^  only  one  saddle  point  and  it  is 
regular.  For  0^  < 0 0^»  the  m^  branch  has  three  saddle  points  labelled  a, 

b,  and  c.  One  or  two  of  these  saddle  points  are  irregular,  because  they 
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lie  on  the  branch  cut  corresponding  to  the  = 0 branch  point  and  the 
approximation  in  equation  (53)  is  not  applicable.  If  < (p  < tx / 2.,  one 
irregular  saddle  point  exists  fcr  the  branch.  In  the  nondispersive  case, 
it  is  known  that  such  irregular  saddle  points  are  associated  with  head  wave- 
fronts  and  two-sided  quasi-equivoluminal  wavefronts.  In  the  sequel,  we  show 
that  analogous  behavior  is  true  in  the  dispersive  case. 

Sketches  of  Cagniard-deHoop  contours  for  0 = 20  and  70  deg  are  shown 
in  Fig.  7.  Im(f.)  is  zero  and  Re(f.)  is  increasing  corresponding  to  increasing 
T / r on  these  contours.  The  Cagniard-deHoop  contours  go  through  the  saddle 
points  and  are  also  paths  of  steepest  descent.  For  a saddle  point  on  a branch 
cut,  the  steepest  descent  approximation  in  equation  (53)  has  to  be  rederived. 

An  irregular  saddle  point  also  occurs  in  the  homogeneous  isotropic  case  [8] 
when  p is  greater  than  the  critical  angle,  which  causes  the  equivo luminal 
wave  to  become  two-sided  behind  the  head  wave.  A two-sided  wave  is  a wave 
that  has  a disturbance  before  and  after  the  arrival  time.  The  present  irregular 
saddle  points  are  treated  in  a manner  similar  to  the  homogeneous  isotropic 
case,  and  the  resulting  changes  are  given  in  Table  I for  the  four  possibilities 
shown  in  Fig.  7 at  points  a and  b for  (p  - 20  deg.  The  coefficient  in  front  of 
the  steepest  descent  approximation  in  equation  (53)  will  change  from  1 to  - 1 
when  the  paths  are  reversed  as  illustrated  by  paths  © and  © . On  paths 
and  , the  coefficient  becomes  -i.  The  other  modification  involves 

whether  the  contour  is  arriving  at  a saddle  point  (paths  © and  0 ) or 
leaving  a saddle  point  (paths  © and  © ).  The  waves  at  these  two  saddle 
points  for  <p  = 20  deg  are  two-sided  because  the  saddle  points  are  on  a branch 
cut.  The  front  portions  result  from  paths  © and  © , whereas  the  back 
portions  result  from  paths  © and  © • This  causes  a sign  change  in  6^, 
where  p is  used  for  front  portions  and  -p  for  back  portions. 

The  contours  in  Fig.  7 are  similar  to  those  in  Kraut's  work  [z].  Cross- 
ing of  the  branch  cut  emanating  from  the  m^  = m^  branch  point  is  done  by 
going  on  to  the  next  Riemann  sheet  of  the  function.  It  is  unnecessary  to  show 
paths  going  up  along  the  cut,  circling  the  branch  point,  and  coming  back  down 
since  the  two  resulting  integrals  are  of  opposite  sign  and  therefore  cancel. 
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A 


Table 


Path 


I*  Coefficients  and  Sign  Changes  in  the  Steepest  Descent 
Approximation  for  the  Qiiasi-Equivoluminal  Two- 
Sided  Wavefronts 


Equation  (5  3) 
coefficient 


Sign  Change  in  6. 


Type  of  wave 


d 


The  six  separate  contributions  to  the  vertical  displacement  v for  t = 5 
and  ~ 45  deg  are  shown  in  Figs.  8-13.  The  nondispersive  solution  (dashed 
curve)  is  included  for  comparison.  These  solutions  would  partially  overlap 
because  the  arrival  times  are  close  together;  however,  they  are  presented 
separately  in  order  to  emphasize  the  difference  in  amplitudes  (note  that  the 
vertical  scale  factors  are  different)  and  to  illustrate  the  effects  of  geometrical 
dispersion  on  the  separate  contributions.  The  effect  of  dispersion  is  par- 
ticularly dramatic  in  Figs.  12  and  13,  where  the  two-sided  wave  that  cor- 
responds to  saddle  point  a is  practically  eliminated.  This  two-sided  wave 
at  point  a is  quasi-equivoluminal  and  is  traveling  in  a direction  that  corres- 
ponds to  or^  - 7 6 deg  which  is  almost  parallel  to  the  laminations.  Previous 
investigations  [4,  5]  revealed  that  these  shear-like  waves  traveling  parallel 
to  the  layering  are  highly  dispersive. 
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Figure  10.  Vertical  Displacement  Due  to  Two-Sided  Quasi- 
Equivoluminal  Wave  Corresponding  to  Saddle 
Point  b;  Back  Part 


d 
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10  v/P 


V.  HEAD  WAVE  CONTRIBUTION 


The  head  wave  is  a reflected  quasi-equivoluminal  wave  generated  by 
the  grazing  incidence  of  the  quas i-dilatational  wave  with  the  free  surface. 

It  consists  of  a planar  wave  that  joins  points  3 and  e in  Fig.  6.  The  head 
wave  is  tangent  to  the  quasi-equivoluminal  wave  at  point  3 and  its  nondis- 
persive  solution  behaves  as  O ( | t - t ' | ^ ' ^)  for  < <f>  < tt/2,  which  is 
rather  weak  compared  with  the  body  wave  singularities.  The  contour  for 
<p  = 70  deg  is  shown  in  Fig.  7;  the  contour  passes  by  the  branch  point  located 
at  Im(q)  = The  portion  of  the  contour  that  is  on  the  branch  cut  just 

below  the  branch  point  corresponds  to  the  region  just  after  the  arrival  of  the 
head  wave. 

In  order  to  estimate  the  effect  of  dispersion  on  the  head  wave,  we 
begin  w'ith 


u(x,  y,  s)/P 


s G ^ A 


exp  (-rsf^)  dq 


(61) 


and  expand  f down  the  imaginary  axis  near  the  branch  point  as  (cf.  Fig.  6.  5, 

[8]) 


+ (q  + f 


= <^ds  ' 


z)/r 


(62) 


= -i/p 


1 1 


where 


^ds/^  = (i/P^i)  sin  0^  + m^ 


COS  0 

c 

q - -i/Pii 


(63) 
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muimimi 


and  since  is  positive  and  increasing  on  a path  of  steepest  descent 


;/r  = 


fz 


p = iq  - 1 /P 


11 


(64) 


q = -i/p 


11 


Angles  with  subscript  c denote  critical  angles  and  correspond  to  point  3 in 
Fig.  6.  Note  that  since 


P/P  ^ = sin 


pm. 


= cos  a 


-i/p 


11 


(65) 


we  can  reduce  t,  to 
ds 


T , = (x  sin  a + y cos  a )/p. 

ds  c c 


(66) 


The  branch  function  is  found  to  have  the  form 


m^  ~ -im^  (p) 


1/2 


(67) 


in  the  fourth  quadrant  just  below  the  branch  point,  where  m^  is  a constant 
and  Gj  is  also  of  thi 
quadrant,  we  obtain 


and  Gj  is  also  of  this  form;  thus,  by  including  the  contribution  from  the  third 


Re/®^r 


Tn , 


tr(x,y,s)/P  = — 


1/  rf: 


•ST 


ds 


q = -i/p 


/■ 


11 


exp  (-sz)  dz 


(68) 
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The  integral  can  be  evaluated  and  the  transform  inverted  to  yield  the  long- 
time result 


Re 

lim  ! sGj^  ^ 

*^1  - ^ds>^^^ 

- -^ds^ 

7T 

s — 0 1 Am^  1 

q = -i/Pii 

(69) 


for  the  nondisper sive  case.  A convolution  then  provides  the  dispersive 
solution  and  the  u displacement  becomes 


u(x,  y,  t)/P  = 


^ Re 

lim 

m^ 

/ a;  \ r^T-v 

3/2 

TTr 

^2 

3/2 

1 P Al  (^c  dp 

CD. 

It 

cr 

(70) 


where  ^'^6  6^  are  defined  in  equations  (59)  and  (60).  In  order  to  obtain 


the  expression  for  v,  G.  is  replaced  by  G.D_.  These  displacements 

-3/Z  ^ i ^ -1/2 

decrease  as  r , whereas  the  body  waves  decrease  as  r . Also,  the 


nondispersive  head  wave  is  smooth  at  its  wavefront,  and  the  convolution 
smooths  it  further.  This  solution  is  not  applicable  when  f^  = 0 at  point  3, 
where  the  head  wave  becomes  tangent  to  the  quasi- equivoluminal  wave. 


VI.  RAYLEIGH  WAVE  CONTRIBUTION 


The  largest  contribution  on  the  surface  in  the  far  field  is  due  to  the 
Rayleigh  wave.  When  y is  zero,  the  Cagniard  contours  encounter  the 
Rayleigh  pole  located  as  shown  in  Fig.  5,  and  by  approximating  the  integrals 
in  the  vicinity  of  the  pole,  it  is  possible  to  determine  the  nature  and  strength 
of  the  Rayleigh  wave  singularity  for  the  nondispersive  case.  A low- 
frequency  expansion  and  a convolution  then  yields  the  dispersive  solution  as 
in  the  previous  sections  for  the  body  waves  and  head  wave. 

The  Rayleigh  denominator  is  first  expanded  as 


The  integrals  are  then  approximated  over  small  portions  of  the  contour  in 
the  fourth  quadrant  along  the  imaginary  q axis  just  above  and  below  the  pole 
The  first  integral  of  equation  (33)  becomes  the  Cauchy  principal  value 


Substitution  of 
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and  completion  of  the  resulting  integral  yields 


I 


I = 

2t7 


P.  V 


/“  -f 

7^ 


-ST 

7f 


dT 


R 


/sG,\ 

U'/ 

q = -i/3j^ 

(74) 


The  long-time  result  is  then  obtained  by  inspection  as 


lim  sG- 


I = 


1 

2tt  s-o  . x/3„)A' 


R' 


(75) 


q = -i/3 


R 


The  integrals  are  also  evaluated  in  the  third  quadrant,  and  the  conjugate 
properties  of  the  integrals  are  used  to  reduce  the  far-field  nondisper sive 
results  to 


” 

- Re 

7T 

lim  ® 

1 

s — 0 A' 

... 

■!  = -'SrO 

= 0 


(76) 


v(x,  0,  t)/P  = — Re 

IT 


lim  ® ^^2^1  ~ ^1^2^ 


s — 0 


q = -i/3 


T - x/3 


RO 


RO 


(77) 


where  the  nondispersive  Rayleigh  wave  velocity  is  given  by 


^RO  ^R 


(78) 
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Note  that  u is  zero  since  the  quantity  in  brackets  in  equation  (76)  is  found 
to  be  imaginary  for  the  parameters  selected  in  this  investigation.  However, 
the  Rayleigh  pole  also  adds  a residue  term  to  the  solution  for  u,  and  this 
can  be  evaluated  by  considering  a small  circular  contour  around  the  pole 
to  yield 


u(x,  0,  t)/P  = Re 


lim  ■ ^2* 

s^O 


q = -i/3 


6(T  - x/3^q) 


RO 


(79) 


An  analytical  expression  for  the  derivative  of  A with  respect  to  q is  given 
in  Appendix  D. 

The  dispersive  solution  corresponding  to  equation  (77)  is  developed 
by  expanding  the  Rayleigh  wave  velocity  as 


1/3^  - ^^^RO  ^R®  ^^RO 


(80) 


and  using  the  convolution  theorem  to  yield 


= i Rp  ! s (G^D^  - ^1^2^  J_  f 

^R  *^0 


v(x,  0,  t)/P  = - Re  . 

IT  I S — ► U 


-1 


A' 


- Ai 


+ p - x/3 


RO 


'R 


dp/ 


Ai 


T - p - x/3 


RO’ 


'R 


(81) 


where 


jl/3 


(82) 
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The  appropriate  dispersive  solution  for  the  residue  term  is  easily  verified 
to  be 


u(x,  0,  t)/P  = Re 


lim  is  (G^  - G^) 


s — 0 


A'X 


R 


q = -i/p 


RO 


. (83) 


Numerical  results  for  the  vertical  displacement  on  the  surface 

[equations  (77)  and  (81)]  are  shown  in  Fig.  14  for  the  same  parameters  as 

the  previous  figures.  The  value  of  in  equation  (80)  for  the  expansion  of 

the  Raylieigh  wave  velocity  was  determined  to  be  approximately  1.03  by 

fitting  a parabola  through  s = 0 and  each  of  the  points  s = 0.  01,  0.  02  and  0.  03. 

This  yielded  three  values  for  Q that  were  then  used  to  find  r at  s = 0 by 

K.  ®r\ 

extrapolation  with  another  parabola.  Note  from  Fig.  14  that  the  dispersion 
eliminates  the  Rayleigh  wave  singularity  and  alters  the  entire  shape  of  the 
response  curve.  Similar  results  are  shown  in  Fig.  15  for  = 0.08  cor- 
responding to  Y = 10. 
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Figure  15.  Vertical  Displacement  Due  to  Rayleigh 
Surface  Wave  (y  = 10) 
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VII.  CONCLUSIONS 


An  approximate  far-field  solution  has  been  presented  for  the  response 
of  a periodically  laminated  composite  half  space  subjected  to  a line  impulse. 
The  governing  equations  were  from  the  effective  stiffness  theory  [4],  The 
solution  consists  of  convolution  integrals  that  arise  when  the  phase  velocities 
are  expanded  for  low  frequency.  As  the  dispersion  caused  by  the  laminations 
is  allowed  to  approach  zero,  the  solution  reduces  to  the  wavefront  expansions 
for  the  corresponding  homogeneous  transversely  isotropic  half  space. 

Numerical  results  were  presented  that  illustrate  the  effect  of  dispersion 
on  the  body  waves  and  the  Rayleigh  surface  wave.  In  each  case,  the  wave- 
front  singularities  are  eliminated  by  the  dispersion,  and  the  displacements 
becoine  continuous. 
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APPENDIX  A.  DERIVATIVES  FOR  WAVE  SURFACE 


The  derivatives  required  to  locate  the  nondispersive  (s  = 0)  wave 
surface  are 


9m 


s = 0 


(84) 


!=0 


(85) 


where  the  primes  represent  partial  differentiation  with  respect  to  q.  The 
derivatives  of  the  B's  for  s = 0 simplify  to 


^5 


s = 0 


s=0 


s=0 


= A^  - ZA^A^/A^^ 


s = 0 


= A3  + A^Ag/Ajj  + AgA^/A^Q 


s =0 


~ A^  - 2AgAg/AjQ 


s = 0 


(86) 


where 


A^^  = 2q[T)Y5^  + (1  - T))  6^] 


A ' = e f 1 
3 m 


^5  = - €^) 


a1  = 2q(T)Y  + 1 - T)) 


Ag  = T1(y  - 1) 


(87) 
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APPENDIX  B.  SECOND  DERIVATIVE  AT  A SADDLE  POINT 


In  order  to  determine  9 ^m/9  at  a saddle  point,  we  begin  with 


*Q(m,  q)  = 0. 


Implicit  differentiation  yields 


and 


0,  q 0,  m , q 


0,  qq  0,  mq  ,q  0,  mm  ,q  0,  m , qq 


However,  at  a saddle  point 


m = -i  tan  <t> 

> q 


leading  to 


, qq  ' 0,  qq  0,mq  ^ 0,  mm  0,  m 


where 


= 2m(2B-B'  - B'B.  - B.B') 
^0,  mq  '3  3 2 4 1 5 


s = 0 


= IZm^B.B.  - 2(B,B.  + B.B.  - B^] 
^0,  mm  14  '24  13  3 


s=0 


(88) 
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(90) 
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The  remaining  derivative  $q  is  found  by  differentiating  the  determinant 


Afm  - 


- A^nni 


A^m  A^ 


s = 0 4 


A^m  - A^ 


which  yields  a sum  of  16  determinants 


0.  qq  (^10^11^ 


4 4 

LE 

s = 0 j = l k=l 


where  ^qjj^  is  with  columns  j and  k replaced  by  their  first  derivatives. 

Note  that  = ^0kj‘  j “ column  j is  replaced  by  its  second  deriva- 

tive. Appendix  A contains  the  first  derivatives  of  the  A's;  the  second  deriva- 


tives are 


A^'  = 4T1y6^  + (1  - Tl)6^] 


= 2(11Y  + 1 - T1) 


J 


APPENDIX  C.  EXPANSION  OF  SLOWNESS  IN  POWERS  OF  s 


The  expansion  for  the  slowness  in  equation  (57)  can  be  obtained  by 
recognizing  that 


m.  = cos  CK:  ^ 

J Uj 


J , SS  I s 

2 I s = 0 


and  since 


”^,ssU  = 0 ■ ^ ,ss ,ml  s = 0 


we  find  that 


C;  - 


-P?  $ 

J_  , SS 


J 2 cos  a-  $ ^'s  = 0 

J „ 

q - q-. 


These  derivatives  can  be  determined  by  evaluating 


$ 


,m 


s = 0 = B,Bj  - b‘)|^ 


*,.sls  = 0 = 'Vl.ss  + - 'BjN.ss  + N®2,s, 


cc  + - 2B,B-  )m‘^  + + B,B_  I 

1 D,ss  5 i,ss  3 3,ss  5 2,ss  2 5,ssis  = 0 


where 
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2 2 

B.  1 « = A.  - 2A,A.  ! ^ A r\  ^ A f\  / 1 n 1 n 

l,sss=0  l,ss  4 4,ss  10  4 10, ss  10’s=0 


i „ = A,  , /A^ . i „ 

2,ss's=0  5 11, ss  ll‘s=0 


®3,  ss  ' s = 0 ^5^9,  ss^^l  1 “ "^5^9 "^1  1,  ss^"^l  1 


-A*AoAj«  /A.^  + AqA.  /Aj«  n 

4 8 10,  ss  10  8 4,  ss  10  s =0 


®4,ss'i  = \,ss  - *ll,ss/^lll=  =" 


2 2 

Br  n “ A C\  A 

5,  SS  s = 0 8 10,  ss  10  s = 0 


''6,ss  = <‘+'l''<'  ■ ^,ss  = ’lA6.ss 

^10,  ss  ° *2^*  ' *1  f , ss  ° 1 
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APPENDIX  D.  RESIDUE  TERM  FOR  RAYLEIGH  POLE 
The  derivative  in  the  denominator  of  equation  (75)  for  small  s is 

A'^gJh^  + g^h'-g'Hj-g^h; 

where 

G]/s  = - iD.  - iqD.'  - m.'  - m'A.  - A^E.'  /s 
J J J J 1 4 j 

D!  = (b'  - 2B  m.m.')/im-  B 
J 1 J J J J 

+ (B  - B.m^)  (im(  B.  + im. 

c ij  Jj 

E'/s  = UA'  D.  + i Ag  D.'  -m'A^)/A,„ 

hVs  = - i e - m.'  D.  A/  - m.  D.'  A.  - A„F.'/s 
J m J J o . J J 6 9 J 

fVs  = (iA  ' - m.'  D.  A - m.  D.'  A_)/A  , 

J 5 J J 9 J J 9 1 1 


4 


1 
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THE  IVAN  A.  GETTING  LABORATORIES 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  tne  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory;  Atmospheric  reactions  and  atmos- 
pheric  optics,  c^iemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser -induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantur  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging:  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory;  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry:  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory;  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  den^ty  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 

THE  AEROSPACE  CORPORATION 
£1  Segundo,  California 
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